Dynamics of 2D pancake vortices in layered superconductors 



A. S. Mel'nikov 

Institute for Physics of Micro structures, Russian Academy of Sciences 
603600, Nizhny Novgorod, 46 Uljanov St., Russia 
melnikov@ipm.sci-nnov.ru 



Abstract 



The dynamics of 2D pancake vortices in Josephson-coupled supercon- 
ducting/normal - metal multilayers is considered within the time-dependent 
Ginzburg-Landau theory. For temperatures close to T c a viscous drag force 
acting on a moving 2D vortex is shown to depend strongly on the conductiv- 
ity of normal metal layers. For a tilted vortex line consisting of 2D vortices 
the equation of viscous motion in the presence of a transport current parallel 
to the layers is obtained. The specific structure of the vortex line core leads 
to a new dynamic behavior and to substantial deviations from the Bardeen- 
Stephen theory. The viscosity coefficient is found to depend essentially on 
the angle 7 between the magnetic field B and the c axis normal to the layers. 
For field orientations close to the layers the nonlinear effects in the vortex 
motion appear even for slowly moving vortex lines (when the in-plane trans- 
port current is much smaller than the Ginzburg-Landau critical current). In 
this nonlinear regime the viscosity coefficient depends logarithmically on the 
vortex velocity V. 

PACS numbers: 74.80.Dm, 74.20.De, 74.60. Ge 
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I. INTRODUCTION 



Recently the peculiarities of magnetic and transport properties of multilayer su- 
perconducting structures attract a great deal of attention. These investigations were 
stimulated, in part, by the discovery of extremely anisotropic high-T c superconductors 
(Bi-Sr-Ca-Cu-O, Tl-Ba-Ca-Cu-O) which may be considered as periodic stacks of two- 
dimensional (2D) superconducting layers. They pertain to a larger class of materials includ- 
ing quasi-2D organic superconductors, chalcogenides and artificial superconducting multi- 
layers Y Ba 2 Cu 3 7 _ s / '(Pr x Y^ x )Ba 2 Cu 3 7 ^, DyBa 2 Cu 3 7 /(Pr x Y^ x )Ba 2 Cu 3 7l Nb/Cu, 
Nb/CuMn, Nb/CuGe, Nb/Al — AlO x A common feature of these systems is the weak 

interlayer Josephson coupling which results in a short effective coherence length £ c for the 
order parameter spatial variation along the c direction (perpendicular to the layers). In 
a broad temperature range the £ c value may be much smaller than the interlayer distance 
D. This fact leads to the quasi-2D character of static and dynamic magnetic properties in 
these compounds. The simplest theoretical model taking account for such a behavior is the 
Lawrence-Doniach (LD) model |§. 

Within this theory one can picture a vortex line passing through a stack of supercon- 
ducting layers (at a certain angle 7 relative to c) as a set of 2D pancake vortices connected 
by Josephson vortices. A single 2D pancake vortex has its singularity (supressed order pa- 
rameter) only in one layer. Note that in a zero external magnetic field 2D pancakes may 
appear in internal superconducting layers due to the thermally activated nucleation of 2D 
vortex-antivortex pairs IPJlQfl. The interaction of 2D pancakes and the vortex lattice struc- 
ture for various magnetic field directions were considered previously in Ref. PHlfiH- This 
model was also used for the description of thermal fluctuations in the mixed state of high-T c 
superconductors [fl7Hl9H • 

The peculiarities of viscous flux flow, flux creep and pinning in anisotropic supercon- 
ductors are also under intensive investigation. In particular, the anisotropy of the flux-flow 
conductivity 07 was investigated previously within the time- dependent Ginzburg-Landau 



(GL) theory with the anisotropic mass tensor [2C-22]. The specific feature of layered com- 
pounds is the additional intrinsic pinning for the vortices which are parallel to the layers. 
The influence of this pinning mechanism on the vortex motion and the peculiarities of the 



flux-flow regime for H close to H C 2 were theoretically studied, for example, in Ref. |2"B| , |2"1 
For weak magnetic fields (H c i < H <C H C 2) the angular dependence of the flux-flow conduc- 



tivity in highly layered superconductors was discussed in Ref. ||22|| . According to Ref. |22 
in a tilted magnetic field in the presence of a current parallel to the layers we have the 
following expression for 07: 

= apyH^ = 0) 
; B I COS7 I v ' 

where o"o is the normal-state conductivity in the plane of the layers. Note that 07 does not 
depend on the orientation of the current with respect to the in-plane magnetic field compo- 
nent. The dissipation is produced by the motion of 2D vortices formed by the component 
B c . The numerical factor fj is determined by the core structure of a single 2D vortex and 
may be calculated using the Bardeen-Stephen theory or within the more exact microscopic 
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models (see, for instance, Ref. [p5||). This simple picture is certainly valid for multilayer 
systems with a negligibly small conductivity 0\ of nonsuperconducting layers (i.e. for su- 
perconductor/insulator multilayers). The problem becomes more complicated if we consider 
the vortex motion in Josephson-coupled superconducting/normal - metal (S/N) multilayers. 
For finite conductivity o\ the spatial distribution of the electric field generated by a moving 
2D pancake changes essentially and the additional dissipation in normal metal layers results 
in the increasing of a viscous drag force acting on a moving 2D vortex. As a consequence 
the viscosity coefficient cannot be evaluated within the simple Bardeen-Stephen model. 

In this paper we investigate the peculiarities of 2D pancake dynamics taking into account 
the nonzero value of the interlayer normal-state conductivity. In Sec. II for zero interlayer 
Josephson coupling we obtain the equation of motion for an isolated 2D vortex in the presence 
of an applied transport current jt. Before discussing the specific features of the 2D vortex 
motion, we will first briefly review some results concerning the motion of a vortex line in a 
3D homogeneous superconductor. The equation for the vortex velocity V in this case has 
the form: 

r)V = x n, (2) 

where n is the unit vector which points in the magnetic field direction (z-direction) and </>o 
is the flux quantum. For T close to T c the quantity rj may be estimated within the simple 
Bardeen-Stephen model or calculated more exactly using the time-dependent GL theory 
(see, for instance, Ref. p5fl). There are two terms in the viscosity coefficient (rj = rji + r/ 2 ) 
accounting for two different dissipation mechanisms: 1) the dissipation due to the relaxation 
of the order parameter magnitude in the vortex core; 2) the dissipation connected with the 
relaxation of the gauge-invariant scalar potential 

TidO ft™„ 
^ + 2^ = ^ VW 
if is the electric potential and 9 - the phase of the complex order parameter which is assumed 
to be the function of the quantity R = r — Vt, where r = (x, y). The expressions for r]i and 
T)2 may be written as follows: 

Vl=Vo jR(^\ dR : = 0.279770 (3) 
o ^ ' 

l - J p 2 (VW) (VW - j?) d 2 R = ~Vo^J p 2 (W6)pd 2 R (4) 

a u4> H c2 

where pexp(i8) is the equilibrium order parameter for a vortex line and R - the distance 
from the vortex line axis. The function ip satisfies the following equation 



r/2 = Vo a = rj 

7TV 
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where £ is the coherence length. The numerical factor u is determined by the pair-breaking 
mechanism. If the lack of the energy gap is connected with the strong electron-phonon relax- 
ation, then it = 5.79 (Ref. p6[). For superconductors with high concentration of magnetic 
impurities we have u = 12 (Ref. fl2~7"| ) and a = 0.159 (see Ref. f28|j29"| ). 

Both the dissipation mechanisms take place for a moving 2D pancake. The first term in 
the viscosity coefficient ( rji) is determined only by the spatial distribution of p(R) in the 
film plane and may be calculated in analogy with the case of a 3D vortex line (Eq. (|3|)). The 
second term (f? 2 ) depends essentially on the /i potential distribution. For a moving 3D vortex 
line \i does not depend on z and decays at a characteristic length ~ £ from the vortex line 
axis. For a 2D pancake moving in the n-th superconducting film of the multilayer structure 
the asymptotic behavior of the /i potential is essentially different. The normal currents in 
nonsuperconducting layers lead to the penetration of the field E generated by a moving 2D 
vortex at a finite length along the z-axis normal to the layers and, as a consequence, to the 
decreasing of ip in the plane z = nD (D is the distance between superconducting layers). It 
will be shown in Sec. II that at large distances R from the 2D vortex center fi(R, z = nD) 
decreases as R^ 1 to the distances 

8(T C - T)e 

ri ~ L„ = 



where £ is the coherence length in superconducting layers. Such an extremely slow decay of 
/i(i?) results in the logarithmic divergence of the integral in the expression (|j) for r] 2 (since 
(VV#) ~ R^ 1 and p{R) — 1 for R 3> £). This divergence is cut off at a large length scale 
L v . 

These arguments show that in a 2D pancake motion equation we should take into account 
not only the terms linear in V but also the terms of the order of Vln(L v /^) = Vln(V c /V), 
where V c = 8(T C — T)C,/irh. From this point of view the situation here is similar to the one 
taking place for a vortex line in a neutral superfluid p0r3TJ. The equation of motion for 



a 2D pancake vortex appears to be essentially nonlinear. The viscosity coefficient depends 
strongly on the o\ value and on the vortex velocity V. These results indicate a breakdown 
of linear response within the employed model. However one may expect linear response to 
hold for more realistic models which consider, for instance, the viscous motion of vortex 
lines consisting of 2D pancakes and take into account the interlayer Josephson coupling. In 
this case the logarithmic divergence of n may be cut off at a characteristic length which is 
either the intervortex spacing or the Josephson length (see Sec. III). 

An isolated 2D vortex is surely energetically unfavourable but the solution of this problem 
is important for understanding of the peculiar dynamics of a tilted vortex line. This dynamics 
is investigated in Sec. III. The viscosity coefficient for a vortex line (and hence the factor rj 
in Eq. (0)) are shown to depend on the angle 7. For field orientations close to the layers the 
qualitative deviations from the Bardeen-Stephen theory become substantial. The nonlinear 
effects in the vortex line motion are also considered. The mechanism of these nonlinear 
phenomena is specific for S/N multilayers and differs essentially from the one proposed in 
Ref. [[32"| , |33f for moving flux lines in 3D homogeneous superconductors. 



In real layered structures the interlayer Josephson coupling may be essential and the 
effect of this coupling on the 2D pancake dynamics should be taken into account. For rather 
small velocities V and large angles 7 the logarithmic divergence of rj discussed above may be 
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cut off at a length scale set by the Josephson critical current J c . The logarithmic dependence 
of TfiV) takes place above a critical velocity determined by the J c value. We will discuss 
these effects qualitatively in Sec. III. 



II. VISCOUS MOTION OF A SINGLE 2D VORTEX 

Let us consider an infinite stack of thin superconducting films of thickness d separated by 
normal metal layers of thickness D ^> d. A good theoretical starting point for the description 
of static magnetic properties of this system is known to be the LD model. For the analysis 
of the dynamic phenomena (for T close to T c ) one may use the simple generalization of this 
model (see, for instance, Ref. |24]]) which can be written in analogy to the time-dependent 
GL model in 3D homogeneous superconductors. It is well-known that such an approach 
corresponds to the real situation for gapless superconductors. 

In this section we will neglect the Josephson coupling of superconducting layers (the effect 
of this coupling will be discussed in Sec. III). For this particular case the general equations 
of the LD model are diagonal in the layer index n. The equations for the order parameter 
p n exp(i8 n ) and the current density in the superconducting layer n have the form: 

~ e^lPn + Pn-Pl- i 2 Pn ( VJn ~ ^ A J ' (5) 



8(T C -T) dt rn ^rny^n ^ 

;r/ ' A + = C 2 V ± (pi ( V ± 9 n - ^A ±l , ) ) (0) 



8(T C - T) rn \ dt h rn J * " V n V he 

he 2 ( 2e \ 

in = Z T?Pn V-iOn ~ ^ A -Ln + cr Ej_ n = j sn + a Ej_ n (7) 



7reA z V nc 

Here (xy) is the film plane , (p n (x, y) = (p(x, y,z — nD), A± n = A±(x, y, z = nD)- the vector 
potential component perpendicular to the z-axis, a Ej_ n - the quasiparticle current density, A- 
the London penetration depth, Vj_ = ■§-), E = — V(f — E^ n = Ej_(x, y, z = nD). 

Note that the superconducting layers are assumed to have a negligible thickness. In 
normal metal layers the expression for the current is: j = a{E. Using the total current 
conservation condition we can obtain the equation for E in nonsuperconducting layers and 
the boundary condition for E z at z = nD: 

divE = (8) 
^j{E z {z = nD + 0)- E z (z = nD - 0)) = — (9) 

In the absence of an applied transport current the phase distribution 9 n for a 2D pancake 
in the film n = is 

v ± e n = ^% , (io) 

r 

where zq, Tq are unit vectors of the cylindrical coordinate system. The magnetic field distri- 
bution for a 2D pancake was found in Ref. JlT],|l2j. In the limit X 2 /(Dd) ^> 1 the expression 
for the vector potential A is: 
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ex P I - I z I \/th I I ~\/ 2 Z A t r I I z o x r o, (11) 



AD V AD 



where A = 2X 2 /d. One can see that A at z = appears to be small in comparison with the 
value if dD/X 2 <C 1. This fact results, in particular, in the logarithmic divergence of a 
2D pancake energy and in the logarithmic interaction of 2D pancakes in the vortex-antivortex 
pair for arbitrary separations [|liyi2"l . Below we will consider only the case dD/X 2 <C 1 
since this condition may be easily fulfilled in real superconducting superlattices at least 
for temperatures close to T c . This approximation allows us to neglect the magnetic field 
generated by the moving 2D vortex in Eqs. ©-(0)- For the analysis of the problem of a 2D 
vortex motion we choose the coordinate system R = r — V£ moving with the vortex velocity 
V and assume that the functions p n ,9 n ,(p depend on r and t in the following way: 

Pn = Pn(r-Vt), <p = <p(r-Vt,z), 9 n = 9 nv (r-Vt)+9 t n , (12) 

where 9 l n is the phase distribution corresponding to the transport current 

The quantity 9 nv meets the conditions: Vj_9 nv — > if | R |— ► oo and 

rotV ± 9 nv = 2vr(5 n o5(r - Vt)z (14) 
The functions p n , 9 nv , <p may be obtained from the following set of equations: 

VV ± p„ = fV\ Pn + Pn -pl- i 2 p n {V ± 9 nv + K n ) 2 (15) 



8(T C - T) 



nh -p 2 n (-VV ± 9 nv + %n)= e 2 V ± [pl (VJ nv + K n )) (16) 



8(T C -T) rn V ^ nv h 
-^'Xz = nD + 0) - <p' z (z = nD- 0)) = -A<p n + ^p 2 n (p n - ^WJ™) , (17) 

he 2 

a = a 1 /a , u 



32AV (T C -T) 

The length £ / y/u is the electric field penetration depth. As was pointed out in Introduction, 
the coefficient u may be calculated within the microscopic theory (see Ref. [26,27|). In 
normal metal layers the equation for the scalar potential <p is: 

Ay? = (18) 

We examine at first the 9 nv and (p distributions for | R £. In this case one can put 
p n ~ 1. The functions a n = V9 nv and p> may be obtained from the system of equations 
(14), (16)- (18). Let us further introduce the 2D Fourier transform: 

<p(k, z) = J ip(R, z)exp(ik ■ R)d 2 R (19) 
a n (k) = J a„(R)exp(ik • R)d 2 R (20) 
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In the Fourier representation we get the equations: 

-^ k " a " = 8(7^(1^ - V - a ") (21) 

-ik x a n = 2n5 n0 z (22) 

-J^-k 2 m = 0, for < z < (n + 1)D (23) 

-(^(z = nD + 0) - ^(z = nD- 0)) = (A; 2 + u/e^n ~ • a n (24) 

For n/} < 2 < (n + 1).D the solution of Eq. ( |23| ) has the form: 

(/9 = c n exp(kz) + d n exp(—kz) 

Using Eq. (|24| ) and the condition </?(| z |— > oo) — > we may calculate the coefficients c n ,d n 
and obtain 

_ 8(r c - t)( 

7rn 

2ey9 n k/^ — 27rV^^ n ok x z — 27ri5 n oY x z 
a " = k • V - i£V c k 2 

The Q value is determined by the equation: 



(26) 



Let us consider the k-domain L~ x <c| k |<C where 

8(T C -T)£ 2 y c e 
irhV V ' 

This k domain correspondes to the distances £ <C| R |<C Lt> and appears to be rather large 
for all velocities V -C V c (we will show below that the last condition is fulfilled for all jt 
values which are much less than the GL critical current j c ). In this domain in the linear 
approximation in V we get the following expressions for if and the phase gradient in the 
layer n = 0: 

ink (V- [z x k]) 

^ ~ ~e k 2 F(k) (28) 

_ 27TZZQ x k 2vrk(V ■ [z x k]) / 1 \ 

a ° " k 2 + Vtfjfc* I f5)J ( } 

F(k) = ((1 + A; 2 ^) 2 + 4A; 2 £>V + 4& J Ds(l + k 2 £ 2 /u)cth(kD)) 1/2 ; s = ^ 
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We consider here only the case of small V values for which the condition 

L v >L m = max[£,4=,AV^>]- ( 30 ) 



is fulfilled. For the domain L m <c| R |<C L v we have: 

»(V.[«.xR.]) (31) 

™ - t + «b f 1 - tct) ( v x z °<" ) - R » (R ° ' |v x (32) 

where C is the constant of the order unity. 

The asymptotic behavior of the n(R, z — 0) potential in the linear approximation in V 
for R 3> L m may be defined from Eqs. (pID,(p2|). At large R values fi(R,z = 0) decreases 
as R~ l to the distances R ~ L v (for a nonzero o\ value). Such an extremely slow decay of 
yu(-R) is caused by the normal currents in nonsuperconducting layers. The field E generated 
by a moving 2D vortex penetrates at a finite length 1/Q along the z-axis which results 
in the decreasing of <p in the plane z = 0. For example, for s ^> 1 (which is possible in 
the vicinity of T c and for not very small a values) and R 3> L m we obtain: 1/Q ~ D^fs , 
(p(z = 0) I^VWo- The specific behavior of (J>(R) results in the logarithmic divergence of 
the viscosity coefficient which is cut off at a large length scale L v (see Eq. (|])). 

To derive a 2D pancake motion equation we use the procedure which is similar to the per- 
turbation method developed in Ref. [p5lj31jl - Let us consider the region R -C min[L v , L t ] (Li- 
the characteristic length scale of the jt variation) and search for the solution of Eqs. (|i~5|) , (|T6|) 
in the form: 

Po = / + g; 9 0v = x + r (33) 

where fexp(ix) is the order parameter for a static 2D pancake, g and r- small corrections 
which are of the order of V and Vln(V c /V). Hereafter we neglect the higher order terms in 
V. The equations for the corrections g and r are 

- ^vv/ = e&g + g - M 2 - ^(V X ) 2 - 2£ 2 fV X (Vr + K ) (34) 

f\2evo/h - VVx) = £KV(/ 2 (Vr + K ) + 2fgV X ) (35) 

Note that for V = 0, K = 0, f = the quantities Xp — (pV)x f P = (pV)/ (p- an arbitrary 
translation vector in the plane (xy)) will solve the equations (|34]) , (|35|) for g and r. We next 
multiply Eq. ([54]) by f p and integrate over the area S of the circle | R | < R\ where R\ meets 
the condition L m <C R\ <C min[L v , L t ]. After simple transformations (using Eq. (|35|) ) one 
finds: 

J f p VVfd 2 R- J f 2 Xp (2e^ /h-VV X )d 2 R = 

s s 

= £V c j> R (Vx P (r + K • R) - x P (Vr + K ))dl, (36) 
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where the circuit L encloses the area S. The quantity t(R = Ri) may be taken from Eq. ( fj2|) 
and the potential ipo - from Eq. (p8f). Evaluating the integrals in Eq. (|36D to the logarithmic 
accuracy we finally obtain: 

Here the constant (3 is of the order unity. The term i]q/3 is connected with the dissipation in 
the domain R < £. Note that we assumed that u £ 1 which is consistent with the microscopic 
theory results for gapless superconductors. It is easy to see from Eq. (^) that the condition 
V c /V > 1 is fulfilled for all j t < j c : 

V c L v CT] V C j c 

— = — k, ~ — 3> 1 

For cti — > we have s — > 0, (3 — > a + 0.279 and Eq. fl3"T|) reduces to the usual equation 

rjV = x z (38) 

c 

with the coefficient n = 771+772 determined by Eqs. d^),(f|). For a nonvanishing con- 
ductivity o\ the main contribution to the viscosity coefficient comes from the region 
max[£,,D] <C R <C Since the parameter s = £ 2 cr/(udD) depends on T as (T c — T) _1 
we come to the conclusion that near T c (if s ^> 1) the 2D pancake motion equation appears 
to be essentially nonlinear. We may neglect the constant (3 and keep only the logarithmically 
large terms in Eq. (|3~7D . With the temperature decreasing the logarithmic term propotional 
to ln(V c /V) will give the dominant contribution only for small velocities V. 

All the results obtained above are relevant only to the case of an isolated 2D pancake. 
In real situation we should consider a set of such 2D vortices which form, for instance, a 
vortex line lattice inclined at a certain angle 7 with respect to the z-axis. It is important 
to note that if the distances between 2D pancakes are less than L v then we can not use 



directly Eq. (37) for each 2D vortex. The thing is that neighbouring 2D pancakes may 
change essentially the \x potential distribution in the vicinity of the given 2D pancake ( for 
R < L v ). Therefore, the dissipation in such a system may depend on the arrangement of 2D 
pancakes. 



III. VISCOUS MOTION OF A VORTEX LINE 

In this section we study the motion of a single tilted vortex line in the presence of an 
applied transport current parallel to the layers. At first let us discuss briefly the case of a 
static vortex line which lies in the (xz) plane, contains the origin and forms the angle 7 with 
the z axis. In the layer n the 2D vortex center is at r n = (nDtg , ~f, 0). Therefore, we have a 
discrete sequence of phase singularities situated in the layers: 

rotV9 n = 2n5(r - r n )z (39) 
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The angle 7 depends on the magnitude and direction of the applied magnetic field 
(H = zqHcos'j + x Hsin i y), which determine the energetically favourable vortex structure. 
In the limit of zero interlayer Josephson coupling the field component parallel to the layers 
is unscreened and penetrates the structure unpertubed, while the perpendicular component 



creates a vortex lattice normal to the layers (the so-called "decomposition property" |fL5|j). 
As a consequence, tilted vortices cannot be formed in a tilted external field if the interlayer 
Josephson critical current is zero. It is, therefore, very important to understand the effect 
of a small Josephson coupling. For a static vortex lattice this problem has been previously 
treated within the LD theory (Ref. []I3| , |r5| , |TE|] ). The tilted vortex lattice was shown to be 



energetically favourable for a wide domain of H and 7. For the field orientation close to the 
layer direction this tilted lattice may transform into the lattice parallel to the layers (the 
orientational lock-in transition ||34|| ) or into a new type of vortex arrangement which consists 



of sets of coexisting parallel and perpendicular vortices |T3 . 

The core structure for a tilted vortex line was examined in detail in Ref. [T^|T^]. It was 
found that this structure depends essentially on the ratio Dtg^/Lj where 



dhc 2 
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is the Josephson length. The case £#7 > Lj/D corresponds to the formation of Josephson-like 
vortex cores connecting 2D pancakes. Hereafter we will consider mainly the other angular 
domain tg^f < Lj/D. In this case the in-plane currents are much larger than the Josephson 
currents. This fact provides a possibility to neglect the Josephson coupling in calculations 
of the vortex viscosity 77. 

As was pointed out in Sec. II, the dissipation produced by moving 2D pancakes may 
depend on their arrangement. The validity of this assumption will be demonstrated in this 
section for the particular case of a tilted flux line. For rather small velocities the logarithmic 
divergency of the viscosity coefficient is cut off at a length scale determined by the size of the 
so-called effective core [Wi (i.e. at the length Dtg'y). One can say that the main contribution 



to the dissipation for magnetic field orientations close to the layers comes exactly from the 
effective core region. In this paper we restrict ourselves to weak magnetic fields and consider 
the case when the effective cores in the tilted vortex lattice do not overlap. Such an approach 
is valid if Dtg^ <C a x where a x is the distance in the x-direction between 2D vortices in the 



same plane but pertaining to neighbouring vortex lines. It was shown in Ref. |15| that this 
condition is fulfilled for 

nD z tg^f 

For a moving vortex line we assume the following dependence of the quantities p n , 9 n , ip 
on the variables r and t: 

Pn = p(R-r n ), (p = <p{R,z), e n = e nv + e t n = e v (R-r n ) + e t n , (40) 

Here R = r — Vt and the quantity 9 l n correspondes to the transport current (0). We will 
examine at first the behavior of the functions 9 V and p in the domain | R — r n £ (p n ~ 1). 
At the same time we consider only the regions | R — r n | <C A which gives us a possibility 
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to neglect the vector potential A. Let us introduce the function a = V6 V . Using the 2D 
Fourier transform 



</?(k, z) = J y?(R, z)exp{ik ■ H)d 2 R 
a(k) = f 3L(R)exp(ik ■ R)d 2 R 



and Eqs. fl)-©,© we obtain: 

l(T - T c ) U 

k 2 (f = 0, for nD < z <(n + l)D (43) 



nh /2e \ 
— i^(k ■ a)exp{ik x nDtg^) = — ^p- ( — <p n — (V ■ a)exp{ik x nDtg^) ) (41) 

-ik x a = 2nz (42) 

dz 2 

fa' z (z = nD + 0)-<p' z (z = nD- 0)) = (k 2 + u/e^n - ^(V ■ s.)exp(ik x nDt gi ) (44) 

After the solution of this system we have the following expressions for the phase gradient a 
and the ip potential in the plane z = 0: 

P(k) = 1 - 1 + — ^- + 2s , ^A ch(kD) - cos{k x Dt 91 )) 
\ u sh{kD) J 

Let us now examine the case of small velocities V for which the condition 

L v > max[£, -4=, D, D^/s, Dtg^A 
\/u 



is fulfilled. To solve Eqs. ([T3|) , ([T0|) for the layer n = we use the perturbation method 
considered in Sec. II. After the similar procedure one obtains the equation (|36|) where the 
radius R\ meets the condition 

maxlC,, — =, -D, -Dv^; £%7] ^ Ri <^ -ky 



Using Eqs. (]4"5|) , (JlTf) we may evaluate the integrals in Eq. (^) taking into account only the 
first order corrections in V. One finally has the equation of the vortex line motion in the 
form (|38|) where the viscosity rj is given by the expression: 



V-Vo 



1 ]/i P(k)d 2 k 
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Here (3q is the constant of the order unity which is determined by the dissipation processes in 
the regions | R — r n |< £. We evaluated the integral in Eq. fl47D to the logarithmic accuracy 
(i.e. for the case when the main contribution to 77 comes from the domain | R — r n |> £) 
and obtained: 

l.D <£,/%7<£,u> 1 



2.D<^t,Dt gi ^t,u 




7) ^ rjo \P + U - -^f^Ts) l < Dt w/0 + InVT+s (49) 

3. D > f,t07 < l,u > 1 

^,o(.g + / n ( (1 + 2 f + / ^ )) (50) 

4. L> > ^,^7 > l,u > 1 



The quantity (3 is of the order unity. These results show that for L v 3> Dtg^y ^> max[D,^] 
the logarithmic divergence of the viscosity rj is cut off at the length Dtgj, i.e. at the distance 
between neighbouring 2D pancakes. This fact results in the angular dependence of the vortex 
line viscosity. 

It is natural to suppose that for the case L v < Dtg^y the main contribution to the dis- 
sipation produced by each 2D pancake moving in the n-th layer comes from the distances 
I R — r n \< L v . As a result the viscosity of a tilted vortex line will depend logarithmically 
on V (see also Sec. II). This supposition is proved by direct calculations. For the domain 
L^ 1 <d A; <gC within the linear approximation in V Eqs. (45),(|46|) may be written as 
follows: 

2mz x k 27rP(A;)k(V ■ [z x k]) 
a(k) * + (52) 

y (M = 0)*(l-i W^t ,Xkl) (53) 

e/r 

Let us now assume that L m C L„ <C Dtg^y. Then for the region L m <C R <C L v one ob- 



tains the expressions for a(R) and <p(R, z — 0) which coincide with Eqs. (j31|) , (p2|) . To derive 
the vortex motion equation we should perform the procedure described in Sec. II choosing the 
Ri value according to the condition L m <C Ri < L 8 . The final equation for V coincides with 
Eq. (|3). The condition L v < Dtgq is fulfilled for vortex velocities V > V* = V c £/(Dtgj) 
which correspond to in-plane transport current densities j t ^> j* where 

Note that the results obtained above are valid only for the case j t <C j c . In the angular 
interval tg^f 3> £/D the current density j* may be much less than j c and there exists a 
current domain j* <C j t <C j c where the vortex line viscosity depends logarithmically on V. 
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Taking into account the expression for the electric field E = x B one obtains 

o- uH c2 r] 

where n is the unit vector parallel to B, H c2 - the upper critical field for H || z . The 
flux-flow conductivity 07 may be written as follows: 

_ cr uH c2 r] 



The quantity rj in Eqs. (pq) , (pq ) is given by the expressions (]48[)-(|5TD if j t < j*. For large 
angles 7 (tg^f ^> £/D) and for current densities jt > j* we have 

where Ej_ is the electric field component perpendicular to the z axis. The current-voltage 
characteristics in this case appear to be nonlinear at least for temperatures close to T c 
(s Z 1). Note that for negligibly small o\ values the expression for 07 fl56| ) coincides with 
Eq. ©. 

As was mentioned above, all these results are valid only in the angular interval 
tgj < Lj/D. In the opposite case Dtgj > Lj the Josephson coupling becomes very essential 
for the calculation of r]. The solution of this problem will be described in detail elsewhere. 
The structure of the effective core in this case becomes more complicated ]nj. We may 



distinguish two contributions to the viscosity coefficient (77 = r] 2 D + Vj) which result from 
the motion of 2D vortices (t]2d) and Josephson-like vortices (7^) connecting 2D pancakes. 
As for the term r\ 2 r> (this contribution will be the leading one at least for \l.{xz)\ it has the 
form 

, 2D * m ^ + m y i + ^ j + ^ - __ j In—^— j (58) 

The logarithmic divergence of rj 2 D is cut off at a length scale Lj if Lj < L v . For this angular 
interval the nonlinear effects in the vortex motion appear for V > V^/Lj. 



IV. CONCLUSIONS 



In summary, we have investigated theoretically some qualitatively new features of 2D 
pancake vortex dynamics in multilayer superconducting structures. The ohmic dissipation 
in nonsuper conducting layers results in the increasing of the viscosity coefficient. In the 
limit of zero interlayer Josephson coupling the viscosity coefficient for an isolated 2D vortex 
appears to be logarithmically divergent. This divergence is cut-off at a length scale set by the 
vortex velocity V. As a consequence the dynamic equation for an isolated 2D vortex appears 
to be essentially nonlinear: the viscosity coefficient depends logarithmically on V. The main 
contribution to ohmic dissipation in this case comes not from the core region (R < £) but 
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from the large- R domain (R < L v ). Therefore, for rather small distances between pancakes 
the dissipation produced by their motion cannot be represented as a sum of contributions 
corresponding to individual 2D vortices. This fact was demonstrated above for the particular 
configuration of pancakes forming a tilted vortex line. 

The equation of viscous motion for a vortex line is obtained and the effects of finite 
Josephson coupling are discussed. The viscosity coefficient is found to depend essentially 
on the distance in the x-direction between pancakes and hence on the tilting angle 7. For 
field orientations close to the layers the nonlinear effects in the vortex motion are shown to 
appear for V > V c ^/min[Lj, Dtg'j]. These deviations from the ohmic behavior occur even 
for slowly moving flux lines when the transport current is much smaller than the Ginzburg- 
Landau critical current. In this nonlinear regime the velocity dependence of 77 coincides with 
that for an isolated 2D pancake. Note that the dynamics of 2D pancakes can be strongly 
affected by the interaction with defects of the crystalline structure, which we have ignored 
in this paper. For 2D vortices localized by linear defects in superconducting layers these 
phenomena were studied previously in Ref. |35| . 
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